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j^ ■ fermions. This extends previous calculations of Golterman and Smit, and Daniel 
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defined in Landau gauge with gauge links excluded. We do this for the simplest 
operators, i.e. those defined on a 2^ hypercube, and for tree level improved 
operators which live on 4^* hypercubes. We also consider gauge invariant operators 
in which the "tadpole" contributions are suppressed by projecting the sums of 
products of gauge links back in to the gauge group. In all cases, we find that 
the variation in the size of the perturbative corrections is smaller than those 
with the gauge invariant unimproved operators. This is most strikingly true for 
the smeared operators. We investigate the efficacy of the mean-field method of 
Lepage and Mackenzie at summing up tadpole contributions. In a companion 
paper we apply these results to four-fermion operators. 
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1 INTRODUCTION 

Lattice calculations of hadronic matrix elements are beginning to have a significant impact 
on phenomenology |jl]]. A necessary step in such calculations is the determination of the 
relationship between the continuum operators whose matrix elements we are interested in, 
and the lattice operators whose matrix elements we can compute. The matrix elements of 
lattice and continuum operators differ due to contributions involving loop momenta close to 
the cut-off, i.e. p ~ n/a, where a is the lattice spacing. As long as a is small enough (1/a 
is 2 — 4 GeV in present calculations), the difference can be calculated using perturbation 
theory. In the present work we extend such calculations to the staggered fermion operators 
which are used in the most accurate numerical simulations 0, ^. 

Perturbative calculations with staggered fermions were pioneered by Sharatchandra, 
Thun and Weisz Q, extended to a variety of fermion bilinears by Golterman and Smit 
[§], and extended again to all bilinears by Daniel and Sheard 0. We refer to the latter 
two papers as GS and DS respectively. The present calculation is an extension of the work 
of DS; we follow their notation, and use and extend the technology they developed. We 
recommend that readers unfamiliar with the subject consult the above mentioned references 
for the technical background. 

A single staggered lattice fermion corresponds to Nf = 4 Dirac fermions in the continuum 
limit, due to the existence of sixteen "doublers" . This has two important consequences. First, 
the continuum theory to which one must match the operators has Nf degenerate flavors for 
each type of staggered fermion. Second, to construct local bihnears having all possible spins 
and flavors one must use lattice operators in which the quark and antiquark fields are spread 
out over a 2^ hypercube [0. Since quark and antiquark fields are, in general, on different 
sites, one must make the operators gauge invariant. The most straightforward method is to 
insert the average of the products of gauge links along the paths of minimal length. These 
were the operators considered in DS, but they have not, in general, been used in numerical 
simulations.^ This is in part because they have large perturbative corrections due to the 
fluctuating gauge links. The corrections can, however, be reduced substantially by summing 
up the tadpole diagrams [0, 0. We discuss this issue in detail in the last section. 



The first simulations attempting to extract matrix elements of such operators used op- 
erators in which the average of the product of gauge links is projected back in to the gauge 



group [12|. The projection reduces the fiuctuations to the level of that of a single link. We 
have calculated the corrections for these operators. This turns out to require only a slight 
modification of the work of DS. 

Most recent calculations (e.g. Ref. 0) use operators in which the lattices are fixed to 
Landau gauge and the gauge links are dropped. This removes the problem of fiuctuating 
gauge links, and also makes the operators much more simple to use in practice. The only 
difficulties are the existence of Gribov copies in Landau gauge, and the fact that, formally, 
one cannot define a transfer matrix when using such operators. One can argue, however, 
that these problems do not have a significant effect on the numerical results, at least in the 



continuum limit |jT3|]. In any case, these are not problems in perturbation theory, in which 



Landau gauge is uniquely defined. The extension of the calculation of DS to these operators 



^Two recent calculations have, however, used these operators ||, |j. 



is presented below. A side benefit of the calculation is that we can check that, for gauge 
invariant operators, the results in Landau gauge agree with those in Feynman gauge. 

A major problem with the matrix elements calculated using Landau gauge operators 
is that there are large 0{a) corrections [|n|, |^.0 With staggered fermions the action is 



accurate to 0{a'^), so to remove such corrections one need only improve the operators. It is 
straightforward to improve the Landau gauge operators at tree level, and a simple method 



for doing this has been proposed in Ref. [|T5|. It requires smearing the quark and antiquark 
fields in such a way that the entire operator gets spread out over a 4*^ hypercube. We have 
calculated the perturbative corrections for these improved operators, and we are also using 
them in numerical simulations 0. 

Many of the phenomenologically interesting matrix elements involve four-fermion oper- 
ators rather than bilinears. Again, numerical simulations have used Landau gauge four- 



fermion operators. For these, one can, as for Wilson fermions |TC|, obtain the perturbative 
corrections from those for bilinears. The main difficulty is that Fierz transformations must 
account for the extra fiavors associated with staggered fermions. We present the results in a 
separate publication [|18]. The calculation for gauge invariant four-fermion operators, which 
requires additional work, was done some time ago for the left-left chiral operators by Sheard 
1^. It has recently been repeated by Ishizuka and Shizawa ||T^. Results for the left-right 



penguin operators are presented in our companion paper ||T8|| . 

The outline of this paper is as follows. In Section ^ we introduce the notation and give 
the Feynman rules. Section ^ gives the details of the computations. We collect the results 
in Section ^, and close in Section |^ with a discussion of the results. 

The results of this work have been presented in preliminary form in Refs. pO| and ||T5 |. 
Various minor errors in these earlier presentations are corrected herein. Ishizuka and Shizawa 
have independently calculated the corrections for the gauge invariant and unsmeared Landau 
gauge operators [0, with results that are in complete agreement with ours. 



2 NOTATION AND DEFINITIONS 

In this section we present the definitions and the Feynman rules necessary for the subsequent 
calculation. A number of the results can be read directly from DS, and we do not repeat 
those here. Other results we have taken from DS but re-expressed so as to highlight the 
similarity of the lattice Feynman rules to those of a continuum theory with four fiavors. 
Finally, the rules relating to the use of Landau gauge and the smeared operators are new. 

2.1 Dirac matrices 

Following DS, we use a gamma matrix basis for both spin and fiavor. To enumerate this 
basis we use "hypercube vectors", four- vectors whose components are or 1. These vectors 
are combined using modulo-2 operations, which, if confusion is likely to arise, we denote 



^Matrix elements of gauge invariant operators also have 0{a) corrections. These appear to be of similar 
size to those of the Landau gauge operators |^ . We focus on Landau gauge operators since they are simpler 
to improve. 



with a subscript 2 (e.g. =2, +2)- A general gamma matrix is labeled by a hypercube vector 

is = i^'i^'ii'i!' . (1) 

We use Euclidean space gamma matrices, which are hermitian, and satisfy {7/^,71/} = 25^,^. 
An alternative basis is built out of the complex conjugate matrices .^^ = 7* and is labeled 
similarly 

^F = ^['^M'^!' ■ (2) 

In the products of such matrices we use the notation 

7575" = 755' , ^S^S' = ^SS' ■ (3) 

Useful notation for hypercube vectors is 

S-F = J2S,F,, S,=2Y.S,. (4) 

In manipulating products of gamma matrices we make frequent use of the following relations 

757^7^ = i-f-^lF = {-f^lF , (5) 

7m75 = l[,^+s] V^iS) , 7]^{S) = (-)S-<M^'' , (6) 

757. = 7[m+5]Cm(^), C,(^) = H^^>^'^ . (7) 

2.2 Continuum bilinears 

We label quark fields in a continuum theory with four degenerate fermions using upper case 
letters, e.g. Qa,a- These have a spinor index (here a) and a flavor index (here a), both 
running from 1 to 4. Color indices play no role in the discussion of this section, and we do 
not show them explicitly. A general bilinear is specified by a spin and a flavor matrix 

O'st = Q.,alf^fQp,b . (8) 

To keep the notation as clear as possible, we always use 7 matrices for spin and ^ matrices 
for flavor. It is convenient to combine spin and flavor matrices in to a single 16 x 16 matrix 

(75 ® ^f) _ _ 

C^5F* = Q.,a(75®eF)""'^'QA. = Q(75®eF)Q. (9) 

The second form is a useful abbreviation in which we treat Q as a 16-component column 
vector. For example, the scalar-isoscalar bilinear Q^ aQa,a is written as (75 ® ^p) with 
S = F = (0000), or equivalently as (1 ® 1). 

We also need two other sets of matrices which are unitarily equivalent to (7s ® ^f)- Both 
sets trade the indices {a, a} for a hypercube vector |^ ^ : 

{1s^^f)ab = iTr[7i757B7^] (10) 

= E ibAriis^^Fr'^'ii^Br-, (n) 

a,a,P,b 



ilS^^F)AB = Eli-r^ ilS®^F)cDli-r- (12) 

CD 



^ 



^ 



■<r 



Figure 1: Notation for Feynman diagrams. 



We refer to the former set as being in the "hypercube" basis, while the latter we say is in the 
"momentum" basis. The "hypercube" basis is convenient for numerical simulations, while 
the "momentum" basis is convenient for perturbative calculations. The multiplication rule 
for each of the representations is as for a direct product, for example 



(7s ® Cf) (7s' ® ^f') = (7s7s' ® CfCf') = hss' ® ^ff') ■ 



(13) 



2.3 Feynman rules 

Since continuum fermions are constructed from the neighborhood of the 16 lattice poles at 
Pfj, = nAfj^, it is convenient to use the decomposition 



Pf^=P'u + T^^t^ 



-7r/2 < p' < 71/2 



(14) 



For physical quarks the "small" component p' is close to zero: if Pphys is the physical mo- 
mentum then p' = apphys- While p' is conserved by the action, the propagator is a matrix 
acting on the indices A. As shown in Refs. p2|, |^, it can be written (using the notation of 
Fig. 0) 



G-\p,-q) = G-\p' + nA,~{q' + nB)) 
= ~5{q-p') 



mil ® 1 



AB 



^Esin(g;.)(7M®l) 



AB 



(15) 
(16) 



Here 5{q' — p') is the periodic delta function, which sets q'^ = p'^ (mod 27r), and there is an 
implicit Kronecker delta for the color indices. 

This way of writing the propagator makes the correspondence with a four-flavor con- 
tinuum theory obvious: the propagators have the same form, except that different (but 
equivalent) bases are used. The only significant change is the appearance of sin q' in Eq. ^ 
rather than q' . This means that lattice and continuum propagators differ by terms of Ola"^). 

When using this propagator it is useful to keep two facts in mind. First, the result 
actually holds for all p' and q', i.e. one can lift the restriction that these momenta lie close 
to the origin. Second, the matrices (7^ '^^)ab simply permute (up to a sign) the indices A 
and B. Thus, for a given value of A, the element of the matrix is non-zero (and equal to 
±1) only for one value of B. 

Using the results of Refs. |^, P] the one gluon vertex can be written 



V^{p' + TiA, ~{q' + nB), -k) = -ig T, 5{q' + k - p') cos(g' + k^/2) (7^ ® 1] 



AB ' 



(17) 



where fi is the gluon polarization direction, c is its color, and the momenta are defined in 
Fig. ^ Again this form is valid for any choices of q' and p', and differs from the continuum 
vertex at 0{a^). 

We need the two gluon vertex only for tadpole diagrams. In these, the total momentum 
of the two gluons entering the vertex is zero, and both the gluons have the same color. Thus 
the full vertex (given in Refs. |^, ^) simplifies to 



V, 



1 ^2 



fiu 



-^29 



{T,f6{q'-p')6,, sin(g;) (7;. ® 1) 



AB ' 



:i^ 



where c is the color of the gluons, fi and z/ their polarizations, and the notation for momenta 
is as for the single gluon vertex. This vertex, which is absent in the continuum, is of 0{a'^), 
because there is, accompanying g"^, an overall factor of a (which we have set to 1) which 
combines with the a from the sing'. 
Finally, the gluon propagator is 



Dik) 



l^u 



5rr' 



E.4sin^(A;./2) 



4sin(/c^/2)sin(A;,/2) 
'' " Ep4sin2(fc,/2) 



(19) 



where c and d are color indices, and /i and v are polarizations. 

Here a = corresponds to Feynman gauge, a = 1 to Landau gauge. The latter is defined 
on the lattice by maximizing X]n,^Tr(f/^(r2)). This may be ambiguous in general, but gives 
a unique result in perturbation theory. 

Ghost propagators are not required in our calculation. 



2.4 Lattice bilinears 

As mentioned in the introduction, we use a variety of lattice transcriptions of the continuum 
operator Og'P^ . These all make use of the "hypercube fields" introduced in Ref. 0. To define 
these we divide the lattice in to 2'^ hypercubes in one of the 16 possible ways. Points on the 
original lattice are then specified by a vector y labeling the hypercubes (with all components 
even), and a hypercube vector C determining the position within the hypercube. The 16 
components of the staggered fermion field x fo^^ ^ given y are now collected in to a single 
hypercube field, which we label following Ref. |^ 

(20) 



x{y)c = \x{y + C) . 



In the continuum limit x{y)c becomes equal to the four-fiavor field Q, when expressed in 
the appropriate basis [^ 

x{y)c — > Q{y)c = C2^crQ(^,b ■ (21) 

Thus the lattice operator 



C>sF{y) = J2 x{y)c {is ® iF)cD x{y)D ■ 

CD 



(22) 



has the same fiavor, spin and normalization as O^p = Q{is'^^f)Q in the continuum limit. 
We often refer to Osf by the abbreviated form xils ® ^f)x- We also use (75 ® ^f) to denote 
OsF, when confusion with the matrix that this notation also defines is unlikely. 



S0F 



V^ \ 




Figure 2: Notation for matrix elements of bilinears. 

2.4.1 Gauge invariant bilinears 

In general, the operators defined by Eq. |2^ are gauge dependent, since they have quark 
and antiquark fields at different sites. We consider various ways of making the operators 
invariant under gauge transformations. First, we use the standard prescription of inserting 
the average of the products of gauge links along the shortest paths connecting the quark and 
antiquark. These operators we refer to as "gauge invariant", or GI for short. 

We need the vertex for an insertion of these operators between an external quark and 
antiquark, and either zero, one or two gluons. The insertion without gluons is (DS Eq. 14) 



r(0) 



M^'i^(p,-g) 



-^" CD 



-iq-D 



CD 



DB p-iq'-D 



(23) 
(24) 



Here we use the momenta as defined in Fig. H, and have set y = for convenience. The 
fact that the operator is at a fixed position means that there is no overall momentum delta 
function. There is an implicit trace over color indices. 

For this vertex to be useful in Feynman diagrams, we must write it in terms of matrices in 
the "momentum" basis, for this is the basis in which the flavor symmetry of the propagators 
is apparent. In fact, for physical momenta (g',p' -^ 0) the insertion becomes (75 ®^f)ab^ 
and is already in the momentum basis. This shows that the operator is correct at tree level 
in the continuum limit. In loop diagrams, however, q' and p' do not vanish, and the operator 
cannot be written in a simple form. To make progress, DS use the following result 



Jk-C 



W^"]Y.EM{k) 



\C-M 



M 



where 



EMik) = l[Ue'''^^/' + 



1 /^-JfcM/2 I (_\M^Jk^/2 



(25) 



(26) 



Combining this with the Eq. ^ one can rewrite the insertion as 

M|?(p' + vrA, -(g' + TTfi)) = (He'iP'.-^'^^A J] Euip') E^{-q') ^ 

\ /i / MN 



MSN 



^mfn)ab ' (27) 



which is the form we use for calculations. Thus the original bilinear with spin S and flavor 
F has, when inserted in loop diagrams, components with all possible spins and flavors. 
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The single gluon vertex (see Fig. ^ for notation) has also been given by DS. We find it 
convenient to slightly rewrite their result as 



(28 



-*-" CD 



where the function /i^^ is related to the function /Z^^,) of DS by 



There are no implicit sums in this formula. This form of the vertex is useful because the 
function /i(^^(fc) is independent of /c^. 

The two gluon vertex is, in general, rather complicated. We need it, however, only for 
tadpole diagrams, and for these there is a considerable simplification. We postpone our 
discussion of the two gluon vertex until we calculate the tadpole diagrams. 

2.4.2 Projected bilinears 

One problem with the operators just discussed is that their matrix elements are suppressed 
due to fluctuations of the intervening gauge links. Much of this suppression comes from 
the tadpole diagrams. The quadratic divergence of these diagrams, when combined with 
powers of lattice spacing coming from expansion of [/^, leave corrections which vanish in the 
continuum limit as powers of g"^ (and not by powers of lattice spacing). A solution to this 
problem was suggested by Sheard |]10| , who summed up the infinite series of tadpole graphs. 



Another solution, proposed by Lepage and Mackenzie |TT|, is to use the mean value of the 
intervening gauge links as a renormalization factor for the bilinears. We have considered still 
another option |jT2|, namely to project the averaged path ordered product of gauge links back 



to the gauge group manifold. We refer to the resulting operators as "projected bilinears". 
As discussed further below, the projection amounts to the replacement 

-^ Y: V eM^g f A.dx'^] -^ exp[tg^ ^ / A^dx'^ + 0{g')] , (31) 

^^path path •' ^^path .p^th ■' 

where the ^^0{g'^)" term makes no contribution in the present calculation. This replacement 
only affects bilinear operators containing more than one gauge link. Moreover, it only alters 
the terms quadratic or higher order in A^, so that only tadpole diagrams are affected. We 
discuss the changes in the vertices when we calculate these diagrams. 

2.4.3 Landau gauge bilinears 

We can avoid making the lattice operators explicitly gauge invariant by choosing a specific 
gauge (modulo ambiguities due to Gribov copies when they exist). The simplest method is 
to fix to Landau gauge^, and then use the operator given in Eq. |2^ without any gauge links. 

^This brings the gauge links as close to identity as possible. 



There is an implicit trace of the color indices of the x ^i^i^d x fields. The resulting operators, 
which we refer to as "Landau gauge" bilinears, have the same zero-gluon vertex as the gauge 
invariant operators (Eq. |2^), but vertices involving gluons vanish. 



2.4.4 Smeared bilinears 

All the bilinears discussed so far differ from the continuum bilinears at 0{a). To see this 
consider the result for the zero-gluon vertex, Eq. ^ The 0{a) terms come from the factors 

of Em and E^ f\. It is easy to see from the definition, Eq. ^ that Em{p') ~ ct'^^' , if 
p' ~ 0{a). Since |Mp only vanishes when M = (0, 0, 0, 0), it follows that all fiavor breaking 
is suppressed by powers of a, as we saw explicitly above. The leading corrections occur for 
\M\ = 1 and |A^| = 0, or for \M\ = and \N\ = 1, and are suppressed by only one power of 
a. 

For operators containing gauge links there are additional 0{a) terms, coming from the fact 
that the gauge fields are at varying positions in the hypercube. For Landau gauge bilinears, 
however, there are no such additional corrections, and thus it is relatively straightforward to 
"improve" the operators so that they have, at tree level, only O(a^) corrections. 

The source of the 0{a) corrections is the non- locality of the staggered operators. The 

quark and antiquark fields have slightly different phases since they are at different points in 

the hypercube. To remove the corrections one must smear the fields so that their average 

position is the center of the hypercube. One method for doing this was suggested in Ref. 

1^ : take the standard bilinear of Eq. ^, replace the quark field according to 



xiy)A ^\J2xiy + 2/i[l - 2A,]U = j-,J2xiy + A + 2/i[l - 2A,]) , (32) 

and perform a similar replacement for the antiquark field. The resulting operators are spread 
out over a 4^ hypercube. The Landau gauge version of these operators, which are the only 
ones that we use, then have only 0{o?) corrections in tree level matrix elements. We refer 
to them as "smeared" operators. 

The Feynman rule for an insertion of a smeared operator has the same form as that 
without smearing, except that the functions Em are changed, i.e. 



Mf^ip' + T^A^-q'-nB) = me^^"^""^^^' 1 E^M (pO Efi-q') hMSN0^MFN)AB , (33) 

where 



MN 



Ef^ik) = EM{k) i V tl^UJ . (34) 

Both Em and E^ share the important symmetry property 

K ^ -K ■ EMJk) ^ i-f^EMik) , El^ik) ^ i-f^Elfik) . (35) 

^It might appear the the first factor on the r.h.s. of Eq. p^ would give 0{a) differences between lattice 
and continuum operators. This is not so. This factor comes from the fact that the lattice operator should 
be associated with a continuum operator located at the center of the hypercube, rather than at one of the 
corners as we have chosen. This difference is irrelevant for the calculations we do here, in which p' = q' 
always. 
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Figure 3: Diagrams contributing to the renormalization of fcrmion bilinears. 

This turns out to guarantee that the operator mixing is the same for smeared and unsmeared 
operators. It is simple to see using Eq. ^that, for k ~ 0(a), the corrections are all of O(a^) 



Elfik) = { 0{a-) 

0(al*-^l') 



\M\ = 
\M\ = 1 
\M\ > 1 



(36) 



which demonstrates that the smeared operators are indeed improved at tree level. 



3 CALCULATIONAL DETAILS 



In this section we present an outline of the calculation of the one-loop lattice corrections to 
the various bihnears defined above. Those interested only in the results can skip to Section ^. 
For the gauge invariant operators, the calculation has been done by DS, and we restate only 
those portions of that work which are needed to extend the calculation to other operators. 
We also show some simplifications in the results of DS that we have found. 

The diagrams which contribute are shown in Fig. 0, where they are labeled following the 
notation of DS. The "X", "Z" and "ZT" diagrams are common to all bilinears, while the 
"Y" and "T" (for tadpole) graphs are absent for Landau gauge operators. For calculations 
involving Landau gauge operators one must use the lattice Landau gauge gluon propagator, 
while the choice of gauge is irrelevant for the gauge invariant operators. 

Our aim is to determine the relationship between continuum and lattice operators. To do 
this, we compare their matrix elements with qq external states. We set the external physical 
momenta to zero, so that the components of the external lattice momenta are either or 
TT. We label the momentum of the outgoing and incoming quarks as p = Ait and q = Bn 



respectively. We set the quark mass to zero, since terms proportional to m are suppressed 
by powers of a. 

At tree level, all lattice versions of the continuum operator Osf have been constructed 
to have the same matrix element, 



M'sF{iat) = {-fs®^F)AB- (37) 

The superscript indicates the number of loops. At one-loop, the matrix elements of the 
various operators differ. Each can be written 

Mlilat) = Cp-f^M^ (38) 

Mi = di Idiv M^iilat) + 4 M]{lat) , (39) 

where the labels i and j each run over the 256 spin-flavor combinations. The coefficients 
dildiv give the divergent parts of the corrections, and are the same for all versions of the 
lattice operator. They are defined precisely below. The coefficients c\, give the remaining 
finite corrections, which do depend on the type of lattice operator. The color factor, Cp = 
{N^ — 1)/2N for the gauge group SU{N), is ubiquitous and so we have factored it out, along 
with g"^ / (IGn^) , to define the reduced matrix element Aii. 

The X- and Z-diagrams are infrared divergent; to regulate them we use a gluon mass A, 
the dependence on which cancels when we compare lattice and continuum matrix elements. 
Following GS, we isolate the divergent parts using the integral^ 



^'''~ ' (4E,^P (4E,^^ + W)- ^^°^ 



where we use the definitions 



" ^^""^ I-n 4^ ' ^^ " sin(0^/2) , c^ = cos(0^/2) . (41) 

The result for this integral is given in GS: 

Idi, = - HX'a^) + Foooo -lE + l, Foooo = 4.36923(1) , 7^ = 0.577216... . (42) 

The coefficients of Idiv are labeled di in Eq. ^, and are proportional to the anomalous 
dimensions of the operators. These corrections do not cause mixing amongst the bilinears, 
as is well known in the continuum. 

All diagrams contribute to the finite coefficients c}j, and it is useful to decompose them 
as 

4 = Xij + 5,, (Y. + Z + ZT + T,) . (43) 

This equation incorporates the result that only the X-diagrams cause mixing. The fact 
that the Y-diagrams do not cause mixing is not apparent from the results presented in 
DS, although it does follow from their numerical results. We present below an explicit 



^DS use a slightly different regularization, in which the gluon mass appears in both propagators. The 
result differs by a constant which cancels in the comparison between lattice and continuum operators. 
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Type 



X 



«j 



K- 



ZT 



T, 



GI (Feynman) 


X(F) 


Y(F) 


Z(F) 


ZT(F) 


T(F) 


GI (Landau) 


X(F) +X(L) 


Y(F)+Y(L) 


Z(F)+Z(L) 


ZT(F)+ZT(L) 


T(F)+T(L) 


Unsmeared L 


X(F) +X(L) 




Z(F)+Z(L) 


ZT(F)+ZT(L) 




Smeared L 


X(F,SM)+X(L,SM) 




Z(F)+Z(L) 


ZT(F)+ZT(L) 




Projected 


X(F) 


Y(F) 


Z(F) 


ZT(F) 


T(F,PR) 



Table 1: Labels for the contributions of the various diagrams to the different types of lattice operator. 
Corrections for gauge invariant (GI) operators were calculated in DS using Feynman gauge (F). Corrections 
for Landau gauge (L), smeared (SM) and projected (PR) operators are calculated below. Blank entries 
indicate that the diagrams do not contribute. 

demonstration of this fact. Furthermore, we show that Yi depends only on the number of 
links in the bilinear, A = 15 — Fp. This is also true of the tadpole contribution, Tj. Wave 
function renormalization is spin-flavor independent, so that Z and ZT do not depend on 
i. The coefficients X, Z, ZT, Y and T depend on which type of lattice operator we use, 
as summarized in Table |l|. DS have calculated these coefficients for the gauge invariant 
operators in Feynman gauge (i.e. the coefficients in the top row of Table |l|); we fill in the 
remainder of the table in the following. 

3.1 X-diagrams (Feynman gauge) 

The calculation of X(F, SM), the Feynman gauge part of the X-diagram contribution for 
smeared operators, follows precisely the same steps as presented in DS for the unsmeared 
operators. There is a divergent part depending only on the spin 



dsF{X) = as , lYl li^lylsluli, = (^s Is 



(44) 



^v 



and a finite part 



-j—rX(F,SM) 
MsF 



/ , X{¥,^y[)£.i'^ 5M,N+p+cr {l^_LpMSN(7^A.® iMFN)J<^Q (45) 

ppaMN 

fipaM 



The coefficients cjj can be read off from this equation. There is mixing between different 
flavors, but the spin remains unchanged. The integrals in this expression are 



X(F,SM) 



fi,pa 

M 



-2 



nSM 



ls,s.BF- ^^(0) E-V 



l^pa^MflB 



(47) 



where 



Sf, = sm(^(/)^j , c^ = cos^cpf,) , 



(E4)-^ 



, B={A 



T.tr 



(4^ 
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M 


X(F,SM)l^i^ 


X(F,SM)^^i 


X(F,SM)]^^2 


X(F,SM)l^i2 


00 


-0.34392 


-0.40583 


0.08410 


0.11908 


10 


-0.06584 


-0.01708 


-0.02214 


-0.00446 


10 


-0.03842 


-0.09339 


-0.02214 


-0.00446 


00 10 


-0.03842 


-0.01981 


0.03534 


0.05691 


000 1 


-0.03842 


-0.01981 


0.03534 


0.00878 


1100 


-0.05870 


-0.03294 


0.08410 


0.11908 


10 10 


-0.05870 


-0.13458 


-0.03323 


-0.00820 


10 1 


-0.05870 


-0.13458 


-0.03323 


-0.05550 


Olio 


-0.08962 


-0.02584 


-0.03323 


-0.00820 


10 1 


-0.08962 


-0.02584 


-0.03323 


-0.05550 


00 11 


-0.08962 


-0.11053 


0.02303 


0.00461 


111 


-0.13870 


-0.26084 


-0.07130 


-0.10743 


10 11 


-0.14808 


-0.05791 


-0.07130 


-0.10743 


110 1 


-0.14808 


-0.17988 


0.03534 


0.00878 


1110 


-0.14808 


-0.17988 


0.03534 


0.05691 


1111 


-0.02866 


-0.01419 


0.02303 


0.00461 



Table 2: X integrals for smeared operators, X{F,SM)'^f'^. 



Results are accurate to ±0.00001. 



The symbols "B" and "F" here stand for the functions coming from boson and fermion 
propagators respectively. The expressions for the unsmeared operators are identical except 
that Eff -^ Em. 

We have calculated these integrals both for smeared and unsmeared operators. The 
results for the latter agree with those given in DS. The former are collected in Table ^. The 
components not shown can be reconstructed using the symmetry of the integral under pa 
interchange, and by permuting the indices. 



3.2 X-diagrams (Landau gauge) 

The contribution of the Landau gauge part of the gluon propagator to the X-diagrams is 
straightforward to evaluate. The reason is that, as in the continuum, the contraction of 
the k^^ in the gluon propagator with the 7^ in the vertex cancels with the adjacent fermion 
propagator, in the limits of vanishing fermion mass and vanishing external momenta. Thus 
the fermion propagators disappear, and the diagram reduces to a tadpole. It turns out to 
be diagonal in spin-flavor space. 

We now demonstrate this in more detail for the unsmeared operators. The substitution 
Em — ^ Ef^ converts the calculation to that for the smeared operators. Using the Feynman 
rules given above, we find the following contribution to the reduced matrix element 






H {lpt^MSNua<^^MFN)AB ^^^'^P^f^^Mi(p)ENi-<l))SuCa{-4:SpS„B) . (49) 
jiupcrMN "^ 

Combining terms using s^ = 2SfjC^ we see that the integral is symmetric under the inter- 
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A X(L)a X(L,SM)a 






0.00000 


4.19152 


1 


3.05826 


4.97458 


2 


4.07769 


5.12452 


3 


4.59516 


5.17390 


4 


4.92625 


5.19848 



Table 3: Landau gauge part of X integrals for unsmeared and smeared operators. The results depend only 
on A = 15 - F|2. They are accurate to ±0.00001. 



changes p ^->- fi and u ^^ a. Thus we can make the replacements (7^p C?) 1) - 


-^ 5^p(l ® 1), and 


(7^0- (g) 1) -^ 5i,o-(l (g) 1). The expression then simplifies to 




-^5F ' = - > : ilMSN ® ^mfn)ab / B'' E m {(j)) E n {- (j)) . 

MN ■^'^ 


(50) 



As promised, the fermion propagators have been canceled. Next we notice that the integral 
vanishes by symmetry unless M = N, and so we find the final expression 

MfJ-^ = -ils®^F)ABj:i~f'^'^^^ I B'EMi<P)EMi-<P) . (51) 

This shows that the correction is diagonal in spin and fiavor. 

The integral is divergent for M = 0, and gives a contribution to di of d{X{L)) = — 1 
independent of spin or fiavor, for both smeared and unsmeared operators. The finite parts 
are _ 

X{L)sF = ^(-)^-(^-^) / UmoB^ - B^EM{(t>)EM{-(t>)) ■ (52) 

It is simple to see from the definition of Eu that the integral depends only on |Mp. It then 
follows that X{V)sF depends only on A = \S — Fp: X{V)sf = -^(-^)a- The numerical 
values for these quantities, for both smeared and unsmeared operators, are given in Table |^. 
For the unsmeared operators, it is possible to derive a simpler form for the result. This 
can be done, for example, by inserting the definition of Em into Eq. ^. The results are 

X(L)a=o = , X(L)a>i = / 5^(1 - n,<AQ) . (53) 



3.3 Y-diagrams (Feynman gauge) 

We review here the calculation of Y-diagrams — those in which a gluon joins an external 
fermion leg to a link in the operator. Clearly these diagrams are absent for Landau gauge 
operators. 

Consider the matrix element of the bilinear (75 ® ^f)- Using the Feynman rules given 
above, standard manipulations yield (for this calculation it is simpler to display the result 
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using the "single-bar" spin-flavor representation) 

CDE t^p 

where the integral is 



+ hi-r-'^ils ® ^F)cnilpP ® ^)DE-A-f-''] , (54) 



l''{p,fi,D-C) = {D-C)JbF spc, h^cD{<P)^M¥<Pp{D-C)^] . (55) 

The first term in Eq. |5^ comes from the gluon attaching to the outgoing quark line, the 
second from attachment to the incoming quark. 

Expanding the exponential gives a sine and a cosine term 

exp[iz0^(D-C)^] = i{D-C),,-s^ + c^ . (56) 

We show below that the cosine term does not contribute, so we keep only the sine term. 
Since the function h^ does not depend on the /i-th component of 0, the integral over 0^ 
vanishes unless p = /i- Thus we can sum over E giving 

Mlf^ = E ^^(75 ® iF)cD^^^ j:iD-C)l JBF si h^cn{<P) . (57) 

CD ^ 4 ^ J^ 

From the definition of h^, it is simple to see that the integral depends on the spin-fiavor of 
the operator only through the combination A = \D — C\'^, the number of gauge links. This 
allows us to sum over C and D, and see explicitly that there is no mixing between different 
spin-fiavors 

A^s^ = (75 ® eF)^B^(^)A • (58) 

The results are conveniently expressed as 

Y{F)a = Y{F)a-i + /a , (59) 



where the integrals are 



h = f BFs\ = f B/4, (60) 

^ / BFslc2 , (61) 

J (p 

: / BFslc2C3 , (62) 

: / BFslc2C3Ci . (63) 



h 
h 
h 



Of course, Y[F)/^=q = 0, for there is no Y-diagram if the operator has no gauge links. We 
tabulate the numerical values for these constants in Table ^. We have checked the results 
numerically against those obtained by appropriately summing the terms in Table 3 of DS.^ 

^The signs of the constants, labeled Yj^j^ by DS, were inadvertently left out from their Table 3; our 
results correct this error. 
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A y(F)A y(L)A 

0.00000 -0.00000 

1 6.11653 -6.11653 

2 7.42832 -8.15537 

3 8.02785 -9.19031 

4 8.41043 -9.85250 



Table 4: Y-integrals in Fcynman gauge, and from the Landau gauge part of the gluon propagators. The 
results depend only on A = jS* — Fp, and are accurate to ±0.00001. 

Now we argue that the cosine term in Eq. ^ gives a null result. This is true separately 
for the diagrams in which the gluon attaches to the incoming and outgoing quarks. In each 
diagram, Sp must pair with one of the s^^^ coming from h^n to produce a non-zero result. 
The resulting contribution to I^{p, fi,D — C) is 

^ {D-C),{D-C)J B{<P)F{<p)slcl (64) 

J (p 

X [1 + C0S{{D-C)^<j)^) + COs{{D~C)r<l)r) + 3 COs{{D -C)^<j)^) C0S{{D-C)r(t>r)] 

with fi,p,a,T all distinct. The integral is just a number, depending on A but independent 
of /i and p. Thus the contribution to /^ is a symmetric function of p and p, but with the 
restriction that p ^ p. Such a function contracted with 7^^ obviously gives zero. Thus the 
off-diagonal part of the Y-diagram vanishes. 

3.4 Y-diagrams (Landau gauge) 

The calculation of the Landau gauge part of the Y-diagrams turns out to be very straightfor- 
ward. The extra term in the Landau gauge gluon propagator again combines with the vertex 
to exactly cancel the fermion propagator. This leaves behind a tadpole graph with only di- 
agonal spin-flavor corrections. The contribution is the same for both possible attachments 
of the gluon. The resulting renormalization is 

Mif = -4 E ^H^-'^hs ® ^Fuii-r'' j:iD-c)i r bx /^^^(</>) . (65) 

CD ti 'f' 

The integral again depends only on the number of links in the operator. A, so that we can 
sum over C and D. Working out the cases individually, the results take the same form as 
the Feynman gauge Y-diagram results, 

r(L)A = r(L)A-i + li , Y{L)^=o = , (66) 

with the integrals given by 

If = -4 / B^l = - f B/A, (67) 
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-'2 



-'3 



-'4 



-4 / Bhlc2 , (68) 

-4 / 52s?C2C3 , (69) 

J(j> 

-4 / Bhlc2CsC, . (70) 



Numerical integration yields the results quoted in Table ^. It is straightforward to show 
that these results are related to those from the X-diagrams in Landau gauge by Y{L)^ = 
-2X(L)a. 

3.5 Self-energy diagrams 

The Z- and ZT-diagrams give rise to wave-function renormalization. The lattice symmetries 
are sufficient to ensure that this renormalization does not break the flavor symmetry [^], so 
that the contribution to matrix elements is independent of spin and flavor. DS calculated 
these diagrams implicitly using Ward Identities. We have checked their result explicitly, and 
extended the calculation to Landau gauge. 

We need to calculate the renormalization diagrams with non-zero external physical mo- 
mentum for the quark, p' 7^ 0, and then pick the term proportional to p' in a Taylor series 
expansion of the resulting expression. (Of course, the term independent of p' vanishes for 
massless quarks.) Considerable algebra leads to the following results for the contribution of 
the Z- diagram in Feynman gauge: 

d{Z{F)) = -1, (71) 



ZiF) 



2BFci{2 - cl){l - 2Fsl) + ^ - B^ 



4 



(72) 



Numerically we find Z{F) = -0.90135(1). 

The Landau gauge Z-diagram leads to the simple results: 



d{Z{L)) = 1, (73) 

Z{L) = — = 3.058263(1) . (74) 



The divergent contribution satisfies two necessary conditions. First, since the total divergent 
contribution from the Landau gauge part of the gluon propagator must vanish, and since only 
X- and Z-diagrams are divergent, it must be that d{Z{L)) + d{X{L)) = 0. Second, as is well 
known in the continuum, the self-energy in Landau gauge is finite, so that d{Z)+d{Z{L)) = 0. 
The tadpole self-energy diagrams yield the finite contributions 

ZT{F) = / I , ZT{L) = - / f , (75) 

J4, 2 J4, Q 

Notice that Z{L) + ZT{L) = 0, so that the the Landau gauge part of the gluon propagator 
has no net effect on the finite part of the self-energy diagrams. The total finite part of 
self-energy correction is thus 

Z{F) + ZT{F) = 11.33170(1) . (76) 
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u. 




Figure 4: 2-link paths illustrating the effect of projection. 

3.6 Tadpole diagrams 

The final part of tlie calculation is tliat of tlie tadpole diagrams. These are obviously absent 
for the Landau gauge operators. For the gauge invariant operators, it is here alone that the 
projection of the sum of products of gauge link matrices back in to the gauge group has an 
effect. 

First let us recall the result if no projection has been done. Then for an operator with 
A links, the tadpoles in Feynman gauge give 



T{F)a 



-A 



B 
~2 



-A X 12.233053(11 



(77) 



To complete our check that Landau and Feynman gauges give the same results we need 
to calculate the Landau gauge contribution to the tadpoles. The diagonal (i.e. [i = v) 
contribution from the Landau gauge part of the gluon propagator is just — ;|T(F)a. In 
addition, the Landau gauge part of the gluon propagator allows off-diagonal contributions 
between links in different directions (i.e. ^ ^ v). Explicit evaluation gives for the total 



T{L)^ = -lT{F)^-4 f B' 



X 






A = 0,1 


■^1-^2 


A = 2 


slsl{2 + C3) 


A = 3 


sfsliS + 2C3 + C3C4) 


A = 4 



(78) 



It is easily verified that T{L)^ = — |F(L)a. Since we also have X(L)a = — |F(L)a, the 
sum of Landau gauge contributions to the X, the Y and the tadpole diagrams gives zero. We 
have already seen that the infinite parts of the Landau gauge contributions cancel, and that 
the finite parts contributing to the self-energy diagrams cancel. Thus we have completed 
our check that the extra part of the Landau gauge gluon propagator gives no contribution 
to gauge invariant operators. This gives us confidence that our calculation of the extra 
contribution to Landau gauge operators, X(L), is correct. 

Now we turn to the operators containing projected gauge links. We only calculate the 
tadpole diagrams for these operators in Feynman gauge. For 0- and 1-link operators projec- 
tion has no effect. We illustrate the calculation for 2-link operators. Without the projection, 
the average of paths connecting the quark and the antiquark is (see Fig. Hf) 



UU1U2 + U3U, 



1 + ^(^1 + ^2 + ^3- 



-Ai)-^{Al + Al + 2A,A2 + Al + Al + 2A3 A4) + • 



(79) 
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A r(F,PR)A T(F,PR)a/T(F)a T(F, PR)A/r(F)i 



1 


-12.23305 


1.00 


1.00 


2 


-14.59650 


0.60 


1.19 


3 


-14.18207 


0.39 


1.16 


4 


-13.05229 


0.27 


1.07 



Table 5: Tadpole integrals in Feynman gauge, with projection. The results depend only on A = jS — F\'^ . 
Results for T(F,PR)a are accurate to ±0.00001. 

where the A's are Hermitian traceless matrices. In Feynman gauge only the diagonal parts 
of the 0{g^) terms contribute, and one gets the result quoted above, with A = 2. The 
projection replaces this matrix with 

2 

1 + 1(^1 + ^2 + ^3 + ^4) - ^{[A,,A2] + [As,A,]) 

- ^(^1 + ^2 + ^3 + ^4) (Ai+A2 + A3 + A4) + --- . (80) 

o 

The 0{g) terms are unchanged, while the 0{g'^) terms are different. For a pair of parallel 
links, say Ui and f/4, the replacement is 

A\ + Al-^\{Al + Al + AiA^ + AiA{) . (81) 

The crucial point is that while A\ gives the integral /^-B, cross terms like ^1^4 give J^Bci 
which is smaller. Note that the 0{g^) commutator terms in the projected link do not 
contribute in Feynman gauge. 

The tadpole contributions with projections are of course diagonal in spin-flavor space, 
and depend only on the number of links in the bilinear. We find 



T(F,PR)i = -\j B (82) 

T(F,PR)2 = -|/s(l + ci) (83) 

T(F,PR)3 = -i /^(5 + 8ci + 5ciC2) (84) 

T(F, PR)4 = -i / ^(1 + 2ci + 2ciC2 + C1C2C3) . (85) 



The numerical values are given in Table |^. We are interested in the ratios by which the 
projection reduces the size of the tadpoles, T(F, PR)a/T'a. We are also interested in the 
size of the projected tadpoles relative to the single link tadpole, T(F, PR)a/T(F)i. Both of 
these quantities are also given in the Table. The results show that the projection is largely 
successful in reducing the size of the tadpole contribution to the level of that of a single link, 
i.e. it reduces the penalty for introducing extra gauge links. 
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4 RESULTS 

To relate lattice and continuum operators, we must also do the one-loop calculation in the 
continuum. Since the continuum bilinears are local, only the X- and Z-diagrams contribute. 
The tree level matrix elements are, by construction, the same as those of the lattice operators 

M'spicont) = {-fs ® ^f) , (86) 

as long as the lattice results (Eq. ^) are expressed in the appropriate basis. The one-loop 
result is 

Mlicont) = Cpj^d'^'^'M^iicont) , (87) 

dr' = icxs-l)i-\n{\yfi') + l)+ts (88) 

where as, defined in Eq. ^, is (4,1,0,1,4) for spin tensors {S,V,T,A,P), X is the gluon 



mass which regulates infrared divergences as in Eq. ^ fx is the renormalization scale, and 
ts depends on the ultraviolet regulator. We present results for four such regulators, all of 
which have been used in previous calculations, and which differ in their treatment of 75: 



A. Dimensional reduction [^ using the "easy subtraction" scheme of Bernard, Draper 
and Soni [El (DREZ); 



B. Dimensional reduction using the subtraction scheme of Altarelli et al. [^ (DRED); 

C Naive dimensional regularization with an anticommuting 75 (NDR); 

D. Dimensional regularization using the 't Hooft-Veltman prescription for 75 p6| (HV). 

All four regulators use modified minimal subtraction to remove the poles at four dimensions. 
The latter three schemes are discussed in detail in Ref. 1271. We find 



( 


0.5, 


0, 0.5, 


0, 


0.5) 


(DREZ) 


( 


0.5, 


0.5, 0.5, 


0.5, 


0.5) 


(DRED) 


( 


-0.5, 


0, 1.5, 


0, 


-0.5) 


(NDR) 


( 


-0.5, 


0, 1.5, 


4, 


7.5) 


(HV) 



^9) 



for spin tensors (S*, V, T, A, P) respectively. 

Combining Eqs. B^ and RSl we obtain the one-loop operator matching relationships 



,2 



or' = Ol"' + Cp^ Y. i^^M Hl^a/n) + dj) Of , (90) 

i 

where the di = {as — 1) are the same as those appearing in Eq. |3^, and 

Cij = 5»j [{c^s - 1) (2 In vr - Fqooo + 1e) + ts] - 4 . (91) 

The dependence on the infrared regulator has canceled, leaving a finite logarithm. We have 
chosen the scale of this log such that it vanishes if yU = vr/a, which is the typical scale of 
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Operator di Name (a) (b) (c) (d) 



(1®1) 


3 


css 


-31.3573 


-31.3573 


-31.3573 


-8.6674 


(1^6) 


3 


CSP 


46.2255 


10.3456 


0.7775 


1.3880 


(i®U 


3 


CSV 


-0.7743 


-0.7743 


-9.9491 


-1.9304 


(1®U) 


3 


CSA 


32.8444 


10.3273 


-0.4221 


0.9857 


(1 ® C/..) 


3 


CST 


18.3027 


8.4331 


-2.8128 


0.2204 


(7m ® 1) 





cvs 


0.0000 


0.0000 


-9.1748 


-10.6286 


(7/. ® 6) 





CVP 


22.5091 


-0.0080 


-10.7574 


-11.1387 


(7m ®U 





cvvo 


-14.7798 


-14.7798 


-14.7798 


-12.6450 


(7m ® e.) 





CvV2 


10.0407 


0.1711 


-11.0748 


-11.1052 


(7m ® Ci^5) 





CVA2 


10.0588 


0.1892 


-11.0567 


-11.3259 


(7m ® Cm5) 





CVAA 


34.5008 


-1.3791 


-10.9472 


-11.0055 


(7m ® ^M^) 





CvTl 


-3.3958 


-3.3958 


-12.5706 


-11.8881 


(7m ® ^^p) 





CVT3 


22.2422 


-0.2749 


-11.0243 


-11.0590 


(7m. ® 1) 


-1 


cts 


8.1761 


-1.6935 


-12.9394 


-14.3310 


(7mi' ® ^m) 


-1 


ctvi 


-2.4711 


-2.4711 


-11.6459 


-14.1209 


(7m. ® 0) 


-1 


CtV3 


19.4637 


-3.0533 


-13.8027 


-14.4338 


(7m. ® ^M.) 


-1 


ctto 


-8.5873 


-8.5873 


-8.5873 


-13.3042 


(7m. ® Cmp) 


-1 


CtT2 


7.8601 


-2.0095 


-13.2554 


-14.3758 


(7m. ® ^M<x) 


-1 


CtTA 


31.2592 


-4.6207 


-14.1888 


-14.4700 



Table 6: Results for the diagonal part of renormalization constants at one-loop, i.e. di and cu using the 
DREZ scheme in the continuum. The components /i, v, p and a are all different. The finite parts ca are 
given for four choices of operators: (a) gauge invariant; (b) projected gauge links; (c) unsmeared Landau 
gauge; and (d) smeared Landau gauge. Those renormalization constants not shown can be obtained from 
these results, as explained in the text. The error in the results is no larger than 0.0001. 
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Operator-i Operator-j Name (a) (b) 



(7m ® 6) 


(7m ® ^m) 


CyvM 


3.04128 


0.12152 


(7m ® ^ms) 


(7m ® 65) 


CVAM 


-0.64608 


0.20344 


(7m ® ^M^s) 


(7m ® ipvb) 


CvTM 


-1.48592 


0.25360 


(7mi' ® ^ms) 


ilf^i^ ® ^ms) 


ctam 


-0.67632 


0.08712 



Table 7: Results for the ofF-diagonal part of renormalization constants at one-loop, Cy, for (a) unsnieared 
and (b) smeared operators. The components /j, v and p are all different, and when p appears it can take 
either of the two allowed values. Mixing coefficients related by axial symmetry are not shown. The error in 
the results does not exceed 0.00001. 

momenta in lattice loop integrals. In this way the coefficients Qj are a direct measure of the 
size of the perturbative corrections. We discuss this further in the next section. 

The diagonal matching coefficients are tabulated in Table ^, for the DREZ scheme.^ We 
give the results for only half of the operators, because the mixing coefficients are the same 
for the two operators (75 ® ^p) and (755 ® ^fs)- Thus, for example, css = cpp.0 The 
names we give the coefficients indicate the spin and flavor of the operator, and, if there is 
an ambiguity, the number of links in operator, i.e. A = jS — Fp. Thus Cyri corresponds to 
the vector current with "tensor" flavor having only a single link. 

The mixing coefficients are given in Table 0. As noted by DS, the 1-loop diagrams do not 
cause mixing to the full extent allowed by lattice symmetries |^. Since only the Feynman 
gauge part of the X-diagrams gives rise to mixing, the off-diagonal coefficients are the same 
for gauge invariant. Landau gauge and projected operators. The only difference is between 
smeared and unsmeared operators. 

To illustrate the use of these tables we give two examples. A lattice operator which, in 
the continuum limit, matches onto the vector current (^(7^ ® C,u)Q (/^ 7^ ^) is 

9' 



where there is no summation over any indices. Our second example requires the use of the 
axial symmetry. The lattice operator matching onto the axial current Q{'y^5 ® ^fiu)Q is 



X(7m5 ® ^M^)X + Cf 



1677^ 



CVT3 X(7m5 ® ^m^)X + CvTM J2 X{lf^5 ® ipv)X 

pi^p,p+v 



(93) 



As stated above, we have checked the results of DS for the gauge invariant operators 
(Tables 4 and 5 of Ref. 0). The relationship between our coefficients and those quoted by 
DSis 

cgs = 4 + bi^ {as - 1) (Foooo - IE) . (94) 



^The results quoted previously, in Refs. fl^, [20], have various errors which are corrected here. 
^°For the operators with an even number of links this is due to the axial symmetry [^ . For operators 
with an odd number of links we do not know of an argument why this equality should hold at higher orders. 
The exception is the conserved vector and partially conserved axial currents (cys and cap), for which the 
renormalization corrections should vanish. 
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Using this, it is simple to compare the diagonal coefficients with those of DS, although the 
off-diagonal coefficients require more work, since DS use a different basis for the operators. 
GS have also done the matching calculation for gauge invariant operators of the form 
(I^^f) for all flavors F . They extract from their calculation finite flavor dependent constant, 
Gp (which has no relationship to as defined above), and a constant from wavefunction 
renormalization, r. One can show analytically that these are related to our constants Qj by 

lQii'^{aF -r) = cfl^ - I = ts=i + 6 In vr - i - ck^f . (95) 

We have checked numerically that these relationships hold. 

5 DISCUSSION 

How large, and how reliable, are the one-loop perturbative corrections that we have cal- 
culated? We want to know the answers to these questions for g^^^^ ^ 1, the value of the 
lattice coupling constant used in present simulations. To obtain answers we must choose 
a value for fia, and a scheme for the coupling constant. We address these two choices in 
turn. Throughout this section we refer to the renormalization coefficients calculated above 
as "Z-factors". 

5.1 Choosing fia and g'^ 

The scalar, tensor and pseudoscalar operators depend upon /la because they have non- 
vanishing anomalous dimensions. When these operators appear in the expression for a 
physical amplitude, however, they will always be multiplied by a coefficient function such 
that the product is independent of ^. In this sense the choice of /i, and thus ^a, is irrelevant. 
In practice, however, the coefficient function is only known to a given order in perturbation 
theory, as is the Z-factor, so there is some dependence on /i, suppressed by powers of g^. The 
same applies for the dependence of coefficient functions and Z-factors on the renormalization 
scheme. It is for this reason that we give the results as a function of fia and for various 
schemes. 

It is useful to explicitly show how this works. We do so for the simplest example of an 
operator which does not mix with others. A physical amplitude is then given by 

Aphys oc C{^x){Ocont)^, = C{fx)Z{fxa){Oiat)a , (96) 

where the first expectation value indicates a continuum matrix element with renormalization 
scale /i, the second a lattice matrix element with lattice spacing a. The product C x Z should 
be independent of /i, and of scheme, as long as one calculates to all orders in perturbation 
theory. We assume that C{mH) is known, in a particular scheme, at a heavy scale itlh- For 
example, if the operator is a scalar, then C{mH) is the running mass at that scale. The 



solution of the renormalization group equation for C is then (in the notation of Ref. ||27|| ) 



Cii^) = Cimn) 



9 imn) 



^ + ^ 16.^ '^ ^Wo-W^""^' 
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(97) 



Here g{fi) is the running coupling constant in the scheme being used, while /3„ and 7*^"^ are 
the n'th order terms in the beta-function and anomalous dimension, respectively. If we use 
another scheme to define the operators, then the form of the result is the same, but C{mii) 
and 7^^-' would differ. Given the result in Eq. ^, one can show that the Z-factor must take 
the following form in order that the product C x Z be independent of /i: 

Z(H = 1 + C^ [-7^°) M/xa/TT) + cl + f ^) [i7^°)(7(°)-2/3o)ln2(/iaA) 



1677^ 



167r^ 



- (7^'^ + (7^°^ -2/5o)c) ln(/ia/7r) + c'] + 0(/) . (98) 

This form applies if the same definition is used for the coupling constant as in Eq. |^. 
Scheme dependence enters through the constants c and c', and through 7'-^''. Our result (Eq. 
|90|) yields the well known results 7*^°^ = — 2Ci?(ii, and also gives the value of c (c = CpCu for 
those operators which do not mix). 

To obtain the physical amplitude we must multiply Eqs. ^ and ^ We then see the well 
known result that the terms multiplying 5'^(/i) in the second parenthesis in the expression for 
C(/x) are of the same order as the 0{g^) corrections to Z{fi) ("1-loop matching requires 2- 
loop anomalous dimensions"), li g'^{mH) << fl'^(/^), so that we can ignore terms proportional 
to g'^{mH), then the full first order correction to the leading logarithms is 



9M 



-7^'^Mn(/ia/7r) + c + 7^"^A/2/3o'-7^'V2/5o • (99) 



This combination must be scheme independent, as long as we always use the same definition 
of the coupling constant, so that the leading logarithmic term in C{fi) is fixed. Thus the 
scheme dependence of c, which enters through the constant tg as shown in Eq. |91], must 
cancel that in 7^^Y2/5o-H This means that, to evaluate the size of the corrections, we need to 
know both c and 7^^-' in a particular scheme, in addition to the scheme independent results 
for /5o, (3i and 7(0). 

Unfortunately, while we know 7*^^-' for scalar and pseudoscalar bilinears (where the three 
loop result is given by |3^), and for vector and axial currents (for which, in NDR, 7 = to 
all orders), we do not know the result for the tensor. For both the scalar and pseudoscalar 
7(0) = -QCf, and 7(1) = -Cp(101 - 10A^//3) in the NDR scheme.0 Setting Nf = 
(appropriate for the quenched approximation), and using the quenched values /3o = H and 
/?i = 102, we find that 7(°)/?i/2/3o — 'y^^^/2f3Q ^ 2.8. For the tensor the first term alone gives 
?a 1.1. These numbers are small compared to the values of c = CpCa in Table P, and we 
choose to ignore them when discussing the size of the corrections. This allows us to treat 
the tensor bilinear, for which we do not know 7^^-', in the same way as the others. We choose 
to use the values of the cu in the DREZ scheme, although the conclusions are essentially 
unchanged in the other schemes. 

^^The scheme dependence in C{mH) will cancel that of the term proportional to g^{mH) in C{ii). As an 
incidental comment we note that our results for ts for the vector and axial currents allow us to deduce the 
two-loop anomalous dimension for the V — A current in the HV scheme. The result agrees with that of Ref. 

^^We are assuming flavor non-singlet operators, which differ from singlet operators at the 2-loop level. 
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As an aside, we note that it follows from the previous discussion that we can deduce the 
value of the 2-loop anomalous dimension on the lattice from that in a continuum scheme: 

7/2 = 7^1 - 2/5oC + 2/507^°) HKomlf^iat) ■ (100) 

This relationship involves the ratio of A-parameters because the coupling constant is different 
in the lattice scheme. In fact, for the operators under consideration here, we could just as 
well dispense with the continuum scheme altogether. This is not true, however, for the main 
application of our results. This is the calculation of Z-factors for the four-fermion operators 
which result when the W and Z are integrated out of the electroweak theory. Since we cannot 
put the weak interactions on the lattice, as the theory is chiral, the coefficients C(m/f ) must 
be computed in the continuum, and we must match the resulting operators with those on 
the lattice. 

Now we return to the issue of which value of /xa to use when estimating the size of the 
corrections. The value should be chosen so that the higher order terms in Z do not contain 
large logarithms proportional to [f7^1n(/ia)]"'. To do this we should match the lattice cut-off 
~ -K /a with the continuum renormalization point /i. Our choice, as already noted above, is 
to use /i = vr/a. A value of, say, /x = l.Tvr/a would be equally reasonable, and would shift 
the coefficients Ca by 21n(1.7)(ii ~ di. This is smaller than most of the coefficients, so the 
ambiguity in /z is only of minor practical significance. It is also important to note that the 
change in cu depends only on the spin of the operator, so that a change in /x does not alter 
the large difference between css and csp in column (a) of Table ^. 

We must also choose a definition of the coupling constant in the expression for Z, and 
the scale at which to evaluate the coupling, which we call /i'. Different choices give dif- 
ferent higher order constants c'. Lepage and Mackenzie argue that higher order constants 



are minimized if one uses a continuum- like scheme such as MS |Tl|]. They also suggest a 
prescription for choosing yu', which would, in general, give a different result for each opera- 
tor. We have chosen to simplify the discussion by choosing the same scale for all operators, 
and taking this to be the same as that in the explicit logarithms: fi' = fi = ir/a. In fact, 
this choice has been made already in the expression for the Z-factor, Eq. ^ The result 
is that, if the bare lattice coupling is g^^^.^ ~ 1, as in present simulations, one should use 
g'^ifi' = Tc/a) ^ 1.8. Taking this value, the factor by which to multiply the numbers in Tables 
I and is 1.8 x Cf/{IQ7i^) ~ 0.015 ^ 1/66. 

5.2 The size of the corrections 

After this long digression we return to the size of the corrections. Multiplying the numbers 
in Tables ^ and ^by 1/66, we see that the corrections for gauge invariant operators (column 
(a)) vary from about — | to greater than |. The only exception is the mixing coefficients, 
which are uniformly small. The large size and range of the corrections suggests that there is 
considerable uncertainty in the application of perturbation theory to such operators. 

A large part of the perturbative corrections, however, comes from fluctuations in the links: 
there is a systematic increase in ca as the number of links increases. These fluctuations are 
reduced by using projected links, as shown in column (b) by the considerable reduction in the 
range of corrections. The Landau gauge bilinears (column (c)) have no links, and also show 
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a reduced range of corrections, roughly from — ^ to 0. This result, together with the fact 
that these operators are simple to implement, makes these operators attractive for numerical 
simulations. 

Finally, the smeared bilinears have the smallest range of corrections, roughly from —0.2 
to 0. There is very little variation between the vector and tensor operators. In addition, 
the mixing coefficients are all smaller than 1%. Thus, to a reasonable approximation, these 
operators are simply slightly rescaled versions of the corresponding continuum operators. 
This makes them good candidates for further numerical study. 

5.3 Summing tadpoles 

Lepage and Mackenzie show that tadpole diagrams are the main source of the large difference 



between the bare lattice coupling and g^ [|rT|. They suggest a mean field method for 
removing the dominant effect of tadpole diagrams. They have applied this to pure gauge 
theory and to Wilson fermions. Here we apply their method to staggered fermions, in order 
to see how well it reduces the fiuctuations in the operators. 

The idea is that the trace of gauge links f/^ fiuctuates around a value uq, rather than 
1, because of the tadpole diagrams. It is the rescaled links U^/uq which one should expand 
around unity, and use in perturbation theory. In particular, the staggered fermion action 

S = ^[mx{n)x{n) + lJ2Vf^in)xin)[U^{n)xin + ^)-U^{n- fiYxin- fj.)]\ , (101) 

n \ ^i ) 

should be rewritten in terms of ?/' = ^Ju{fx 
S = Y^ (^{m/uo)ilj{n)ilj{n) + 

n 

lY.Vt.{nmnW,{n)/uomn + /x) - {U,{n - i^)yuo)ij{n - /x)] ) . (102) 
^J. ) 

The claim is that V^ is better matched to the continuum quark field than is Xi ^i-nd that one 
should use tl) to construct lattice operators. Similarly, if such operators contain links, then 
they should be the rescaled links U^/uq. For the bilinears we are studying, the prescription 
is then to multiply the naive lattice operator by Uq~'^^ , where nu is the number of gauge 
links. For the gauge invariant operators riu = A, for the Landau gauge operators riu = 0, 
while for operators with projected gauge links we assume that riu = 1 — 5a,o- 

In addition, we should subtract from the perturbative corrections Cij the parts coming 
from the tadpole diagrams. This is because such diagrams are the lowest order contribution to 
Uo, and thus have already been taken into account. The actual value of mq is to be determined 
non-perturbatively, in effect summing up tadpoles (and possibly some other effects) to all 
orders. The precise amount to subtract from Cij is weakly dependent on the definition of 
Uq used. We opt for the definition uq = {^Ti{U)l), the trace of the link in Landau gauge. 
In perturbation theory this means that mq = 1 — CF(5'^/167r^)9. 17479 + 0{g'^). This choice 
has the advantage that it amounts to dropping from cu the contribution from all tadpole 
diagrams (ZT and T) evaluated in Landau gauge. 
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Name (a) (b) (c) (d) 



css 


-22.1825 


-22.1825 


-22.1825 


0.5074 


CSP 


18.7011 


10.3456 


9.9523 


10.5628 


CSV 


-0.7743 


-0.7743 


-0.7743 


7.2443 


CSA 


14.4948 


10.3273 


8.7527 


10.1605 


CST 


9.1279 


8.4331 


6.3620 


9.3952 


cvs 


0.0000 


0.0000 


0.0000 


-1.4538 


CVP 


4.1595 


-0.0080 


-1.5826 


-1.9639 


CVVQ 


-5.6050 


-5.6050 


-5.6050 


-3.4702 


CvV2 


0.8659 


0.1711 


-1.9000 


-1.9304 


CVA2 


0.8840 


0.1892 


-1.8819 


-2.1511 


CVA4 


6.9764 


-1.3791 


-1.7724 


-1.8307 


CVTI 


-3.3958 


-3.3958 


-3.3958 


-2.7134 


CVT3 


3.8926 


-0.2749 


-1.8495 


-1.8842 


cts 


-0.9987 


-1.6935 


-3.7646 


-5.1562 


ctvi 


-2.4711 


-2.4711 


-2.4711 


-4.9461 


CtV3 


1.1142 


-3.0533 


-4.6279 


-5.2590 


CtTO 


0.5875 


0.5875 


0.5875 


-4.1294 


CtT2 


-1.3147 


-2.0095 


-4.0806 


-5.2010 


ctta 


3.7349 


-4.6207 


-5.0140 


-5.2952 



Table 8: Results for the diagonal part of renormalization constants at one-loop, cu, with Landau gauge 
tadpoles removed, as appropriate in the mean field prescription of Lepage and Mackenzie. The notation is 
as in Table 6. 

The same prescription applied to the gauge action itself predicts that the correct coupling 



constant to use in perturbative formulae is f^bare^o illi- For g^^^^ = 1, numerical simulations 



yield (| Tt{U)l) ~ 0.86. Taking this as our estimate of mq, we find that the effective coupling 
should be l.Sg"^, which is the value chosen above. 

To illustrate the effect of this prescription we show the shifted diagonal corrections in 
Table §. Notice that the prescription does not affect the conserved vector current at all, 
since ZT = — Ta=i. It does affect the Landau gauge version of this current, shifting its 
correction to zero. The result of the shifts is, as expected, most significant for the gauge 
invariant operators. The corrections now range roughly from — | to |, possibly small enough 
that perturbation theory is reliable. The range of the corrections is also reduced for the 
projected operators, for which the corrections remain slightly smaller than those for the gauge 
invariant operators. The tadpole subtraction simply shifts the corrections for the Landau 
gauge operators, with the result that they are more evenly distributed around zero. The 
smeared operators still have the smallest corrections. Finally, we note that the prescription 
has no effect on the off-diagonal corrections, because tadpoles contribute only to diagonal 
renormalizat ion . 
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